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abstract 

(^ I We investigate the existence and stability of gap solitons in a double-core optical fiber, 

m 

O ' where one core has the Kerr nonlinearity and the other one is linear, with the Bragg 
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O ■ grating (BG) written on the nonlinear core, while the linear one may or may not have 
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O I BG. The model considerably extends the previously studied families of BG solitons. For 
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zero-velocity solitons, we find exact solutions in a limiting case when the group-velocity 
terms are absent in the equation for the linear core. In the general case, solitons are 
found numerically. Stability borders for the solitons are found in terms of an internal 
parameter of the soliton family. Depending on frequency u, the solitons may remain 
stable for large values of the group velocity in the linear core. Stable moving solitons are 
found too. They are produced by interaction of initially separated solitons, which shows 
a considerable spontaneous symmetry breaking in the case when the solitons attract each 
other. 



INTRODUCTION AND FORMULATION OF THE MODEL 

It is well known that the combination of the Kerr nonlinearity with a strong effective 
dispersion induced by the resonant reflection of light on the Bragg grating (BG) gives 
rise to a vast family of gap solitons, frequently called BG solitons |1|] (in this work, we 
use the term "soliton" in the loose sense, without implying integrability of the model 
where it appears; in particular, it will be shown that interactions between "solitons" 
in a model to be introduced below may be essentially inelastic). A generally accepted 
mathematical model of the nonlinear fiber equipped with BG is the so-called generalized 
massive Thirring model (GMTM) [^. Thorough theoretical investigation of BG solitons, 
an important step in which was the discovery of a class of exact single-soliton solutions to 
GMTM 10], was followed by observation of BG solitons created by a very strong laser pulse 
launched into in a short segment (~ 6 cm) of a nonlinear optical fiber with the resonant BG 
written on it Q. Experimental studies of BG solitons were further developed (including, 
in particular, formation of multiple BG solitons) in Refs. 0]. 

Observation of solitons in such a short fiber paves the way for many potential appli- 
cations, as well as for new experiments aimed at the study of fundamental properties 
of optical solitons. This also makes it relevant to consider more sophisticated nonlin- 
ear systems based on fiber gratings, where the properties of solitons might still be more 
promising. In particular, one can look for solitons in a dual-core system with linear cou- 
pling between the cores, BG being written on both cores or a single one. The case of 
two identical BG-carrying cores was considered in Ref. |]^, where it was found that the 
model gave rise to a bifurcation at a critical value of the soliton's energy. The bifurcation 
destabilizes a symmetric two-component solution, simultaneously generating a nontrivial 
asymmetric soliton. A dual-fiber system with unlike cores is easier to fabricate and may 
offer new possibilities. One of the most interesting dual systems with different cores is a 
semilinear one, where one core is linear. Semilinear dual-core models without BG were 
introduced earlier; both continuous-wave and soliton states in them have been studied in 
various contexts^, 0. 

The objective of this work is to introduce a semilinear dual-core model in which BG 
is written either on the nonlinear core only or on both cores, and to search for solitons 
in it (which makes it necessary, first of all, to explore the system's linear spectrum). 
Following the derivation of GMTM |1[ and of the standard equations for a dual-core fiber 



(see, e.g., Ref. |^) from Maxwell's equations, a general model for the semilinear dual-core 
BG-equipped system can be written as the following set of normalized equations: 

iut + iu^ + \\v\^ + (l/2)|Mpl U + V + ncj) = Q, (1) 



ivt - ivx + 



u\'^ + {l/2)\v\'^ v + u + kiIj = 0, (2) 



i(j)t + ic(j)x + KM + (A + 2/i) i' = 0, (3) 

iipt — icipx + Kt> + (A — in) = (4) 

Here, u and v represent the forward- and backward-propagating waves in the nonlinear 
core, and ip are their counterparts in the linear one, k is the coefficient of linear coupling 
between the cores, while A and /x are the real and imaginary parts of the BG-coupling 
coefficient in the linear core (which is, generally, complex if its counterpart in the nonlinear 
core is normalized to be 1, as is the case in Eqs. (|I|) and (0)). Lastly, the group velocity 
in the nonlinear core is set equal to 1, and c is the relative group velocity in the linear 
core. 

The simplest case giving rise to a novel system is A = yU = (corresponding to the 
linear core without BG), while cross-core coupling k is nonzero. Below, we will always 
set /i = 0; in most cases, A will also be zero, but effects of A 7^ on the solitons' stability 
will be investigated too. Note that, although the present model finds its most natural 
formulation in the temporal domain, it can also be readily interpreted in terms of the 
spatial- domain evolution of the fields in a two-core planar waveguide, BG being realized 
as a system of parallel cores written on the waveguide(s) p. 

It may also be quite interesting to consider a system where the Kerr nonlinearity and 
BG are separated, i.e., with the grating written only on the linear core. The corresponding 
model is obtained from the above equations, dropping the linear terms v and u in Eqs. 
(|^) and (0) and setting A = 1 and /i = in Eqs. (^ and (^). This model, which also 
seems quite promising, will be considered elsewhere. 

Before looking for solitons, it is necessary to analyze the spectrum of the linearized 
system, in order to identify a spectral gap in which BG solitons may reside ]l|]. For a 
linear wave ~ exp {ikx — iut) and setting /i = 0, a dispersion equation for u^k) can be 
obtained: 



a;^- 



1 + 2«:2 + a2 + (1 + c") k^] ^2 + (a - K^y + (c2 - 2CK^ + A^) k^ + c'k^ = 0. 
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Analyzing this equation, it is easy to conclude that the gap does not exist in the present 
model if A < K^ and c^ - A + A^ < (2c - 1)k'^, or if (1 + 2c)'^ {c + c'^ + A^) < k^ < X. In 
all the other cases, a finite gap is present, and BG solitons may exist. In the particular 
case when the linear properties of the two cores are identical, i.e., c = 1 and A = 1, which 
physically corresponds to identical BGs written on them, the gap existence condition 
takes a very simple form, k^ < 1 [|5[. 

A remarkable property of the above-mentioned GMTM equations, to which Eqs. (1) 
- (4) reduce if the additional core is dropped, is the availability of exact single-soliton 
solutions, both quiescent and moving with an arbitrary velocity v, limited by |f | < 1, 
despite the fact that the model is not integrable (except for the unphysical case when 
the self-phase modulation terms are omitted) 0. Here, we aim to find soliton solutions 
to the full system (1) - (4) and investigate their stability and interactions. Solitons with 
zero velocity will be studied in detail, and moving solitons will be presented too. In 
fact, the existence of the solitons with zero velocity (which have not yet been observed 
experimentally in the single-core fiber gratings) is the most intriguing possibility, as this 
implies a possibility of "full stoppage of light" through its dynamical trapping, which is 
especially interesting in view of the recent discovery of "ultraslow light" in ultracold gases 

As for the physical parameters of the system and its soliton solutions, a crucial factor 
is the ratio of the length -Zcoupi of the coupling between the cores and a characteristic 
propagation distance (soliton's dispersion length) Zsoi necessary for the formation of a 
soliton in a single-core fiber with BG. As is well known, the former length in available 
dual-core fibers is, normally, ~ 1 cm, and, according to the experimental data ^, ^ , Zgoi 
is on the same order of magnitude (it is so short, despite the fact that the solitons are 
relatively broad in the temporal domain, because BG gives rise to an extremely strong 
effective dispersion). This circumstance, -Zcoupi ~ -^soi, is quite favorable, as it suggests 
that the interplay between the resonant light reflection on BG, Kerr nonlinearity, and 
linear coupling between the cores may give rise to solitons with fairly novel properties, 
in comparison with both the usual (single-core) BG solitons [0 and solitons in dual-core 
fibers without BG 0. 

In line with the above arguments, other basic characteristics of the new solitons are 
expected to be on the same order of magnitude as those for the recently observed BG 
solitons in the single-core fiber. In particular, the soliton can be generated by a laser pulse 



of the duration ~ 100 ps, having a fairly high peak power ~ 5 W (which is, however, 
still sufficiently far from the optical-breakdown threshold in silica glass) and, accordingly, 
the energy ~ 500 pJ . The soliton to be created will keep, essentially, all this energy, 
self-compressing to the temporal width ^ 50 ps [^ Q. 

Another crucial ingredient of the possible experiment is the necessary length of the BG- 
equipped dual-core fiber. As it was mentioned above, for the successful generation and 
detection of the gap soliton in the single-core BG fiber, a 6 cm fiber was sufficient. In fact, 
the present-day techniques make it quite easy to fabricate a homogeneous dual-core fiber 
of the length ~ 1 m, as well as to write a uniform BG on it. Therefore, the experiment 
may be quite feasible in the fiber whose length is on the order of 100 characteristic soliton 
and coupling lengths (both being ~ 1 cm, see above), which will be more than enough 
for the most precise experiments. 

Thus, experimental generation of the new solitons, to be theoretically studied in the 
present work, is not going to be much harder than the recent experiments reported in 
Refs. and 0]. The only essentially new issue in the experiment may be a question if 
the input laser pulse may be focussed, as usual, on the entrance face of one core only, or 
it is necessary to split it, in a special fashion, between the two cores. Although it may 
be premature here to discuss experimental technicalities in such a detail, we notice that 
having the fiber length much longer than ^coupi ~ 1 cm, see above, will provide enough 
room for the proper redistribution of the power between the cores, so that the experiment 
will scarcely be critically sensitive to details of launching the input pulse. 

The rest of the paper is organized as follows. In section 2 we display exact analytical 
soliton solutions that can be found in the present model with c = 0, and results of 
simulations of their stability, which show that they are stable in a broad parametric 
region. In the case c 7^ 0, soliton solutions can only be found numerically, which is done 
in section 3, together with systematic simulations of their stability. It is found that, 
depending on the value of frequency oj the solitons may remain stable up to a large value 
c = Cmax- At c > Cmax, the solitou bccomcs unstable. This instability, however, does 
not destroy it, but, after shedding some radiation, it evolves into another member of the 
soliton family. 

In section 4 we directly simulate interactions between two solitons placed initially at 
some distance from each other. It is found that the result of the interaction strongly 
depends on the relative phase of the two solitons. In particular, the interaction can easily 



generate moving solitons and leads to spontaneous symmetry breaking. 



EXACT SOLITON SOLUTIONS AND THEIR STABILITY 



Exact zero- velocity soliton solutions to Eqs. (HD - (@) can only be found in the particular 
case c = 0. Starting with the usual ansatz, 



u = U{x) exp{—iu!t), V = V{x)exp{—iut), 
(j) = $(x) exp(— icut), ip = '^{x)exp{—iujt) , 

and following the pattern of the exact GMTM solutions [0] , we find 



(5) 
(6) 



U{x) = 
V{x) = 



(^2_^2_^2^_^^2)2^^2^4nV4 



(cu2 - A2 - n^y 



2^2 I ,,2^4 



(^2 _ ^2 _ ^2 ^ ^^2y ^ ^2^. 



(Cj2 _ A2 - i^^Y 

where 6 = tan~^ (/t^Ai/ {^"^ — A^ — /i^ + Xk^)), and 



1/4 



e+''/'A{x), 



-'^I^B{x) 



(7) 



A (x) = Y^ 2/3 (sin ^) sech {r]x ■ sin - i^/2) , 
5(x) = -./273(sine) sech (r^a;- sin ^ + i^/2) 



(8) 



<^{x) 
^ix) 



KUO 



u+ ^^^:t'^Kv. 



LU"^ — X"^ — /i^ tU^ — A^ — yU^ 



K (A — ZyU) 



:f/ 



KUJ 



:V. 



(9) 



cu^ — A^ — /i^ iu"^ — X"^ — iJ? 
Here ^, which takes values between and vr, is an arbitrary parameter of the soliton 
family. The frequency uo and inverse width r] of the soliton are determined, in terms of 9, 
by equations 

U (u;2 _ A2 - ;,2 _ ^2) 



y^(cu2 - A2 - /i2 + Afi:2) V yU^K 



2^4 



/ 2 \2 2 

sgn [uj — A — /i 



cos 6* , 



(^2 _ ^2 _ ^2 ^ ^^2)^ ^ ^2^4 



nV4 



(10) 



11) 



(cj2 _ A2 - ijfly 

It is relevant to note that these exact solutions resemble those found earlier in a linearly 
coupled system of cubic and linear Ginzburg-Landau (GL) equations |[T^; however, the 
exact solutions to the GL equations exist as isolated ones, rather than in families, i.e., 
they do not contain any arbitrary parameter. 
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Before proceeding to numerical search for solitons in the case c 7^ 0, it is necessary to 
address the stabihty of the exact analytical solutions obtained above. The first nonrig- 
orous stability analysis of GMTM solitons was done using the variational approximation 
TT[ . It was predicted that instability might occur when an internal parameter of GMTM 



solitons 9, similar to that introduced above in Eqs. (^, exceeds a certain critical value, 
which was close to 7r/2. Then, a rigorous treatment of the stability problem for the 
GMTM system, based on the consideration of its linearized version, was developed in 
Refs. [0. It was demonstrated there that the solitons with 9 exceeding a critical value, 
which is slightly larger than 7r/2, are unstable indeed. However, the instability is weak, 
therefore it was hard to observe it in direct simulations. 

In this connection, it should be noted that while the results for the solitons' stability 
in various models, obtained from the solution of the corresponding eigenvalue problem 
for the linearized equations are more rigorous (and usually are technically more difficult) 
than those produced by direct simulations of the nonlinear equations, the latter results 
may be more appropriate for the physical applications. Indeed, if the soliton is, rigorously 
speaking, unstable but the instability is weak (as is the case for GMTM), it may happen 
that neither direct simulations performed for a limited evolution time (or propagation dis- 
tance, depending on the particular system) nor a real experiment in a finite-size sample 
will demonstrate the instability, so that, in terms of real physics, the soliton should be 
regarded as a stable object, in accord with the prediction of the direct simulations, and 
despite the contradiction with the rigorous results. Solitons in the BG fiber may provide 
an example of this situation. In this case, experimental results |P, Q, while being in good 
agreement with direct simulations, have not been able to demonstrate the sophisticated 



instability predicted on the basis of the linearized equations in Refs.|]12[. On the other 
hand, it is necessary to mention that, although the physical value of the soliton's peak 
power in these experiments was quite high, the actually observed BG solitons may still 
be low-intensity ones from the viewpoint of the corresponding theoretical model. How- 
ever, the above-mentioned "sophisticated instability" takes place only for high-intensity 
solitons. Another example that could be cited regarding the fact that solitons' instability 
may sometimes be formal is provided by "spinning" (2+l)-dimensional solitons in me- 
dia with the cubic-quintic nonlinearity. As the analysis of the corresponding linearized 
problem shows, the solitons with the "spin" s = 1 are, strictly speaking, always unstable 
against infinitesimal azimuthal perturbations which destroy the cylindrical symmetry of 



the solitons. Nevertheless, if the size of the "spinning" sohton is large enough, the insta- 
bility may be so weak that the soliton may persist as a fairly robust object over several 
diffraction lengths ^^, thus having a fairly good chance to be observed in an experiment. 



To test the stability of the exact soliton solutions given by Eqs. (^) through (0), we 
simulated their evolution by means of the split-step Fourier algorithm, imposing various 
asymmetric (sometimes nonsmall) initial perturbations. A typical case is displayed in 
Fig. 1, showing that after shedding off some radiation, the perturbed pulse readily evolves 
into a member of the soliton family (in fact, the final soliton in Fig. 1 acquires a very 
small velocity, because the asymmetric perturbation has "pushed" it; moving solitons will 
be specially considered below). In particular, an important finding is that, when A 7^ 0, 
Eq. ([10|) gives rise to three distinct roots for uj, of which only the one with largest |a;| is 
found to produce a stable soliton. On the other hand, there are two different roots for uj 
at A = 0, both leading to stable solitons. 

We have also found that, similar to the GMTM solitons, a fundamental property of 
the soliton family in our extended model is that a stable part of the family is limited, 6 < 
6'max, where 6'maxdepends on k and A. To analyze this in detail, we set A = 0, focusing on 
the simplest and most fundamental case when BG is present in the nonlinear core only, 
and the stability is solely controlled by the coefficient of the linear coupling between the 
nonlinear and linear cores. The stability border inside the soliton family, 6niax{i^), was 
then sought for gradually increasing at a fixed value of k. We started from 6 = vr/12, 
where the exact soliton is definitely stable, until we hit a value 6'max that gave rise to 
instability. The instability, when it sets in, causes straightforward decay of the soliton 
into radiation. We have thus found that 6'max( ^ = 0.01) = 7r/1.7, 6'inax( n = 1) = 7r/2.0, 
and 6ma.x{i^ = 100) = vr/l.S, i.e., the dependence of the stability limit on k is fairly weak, 
6'max being close to that in the single-core model, although the shapes of the exact solitons 
may be really different. 

SOLITONS IN THE MODEL WITH c / 

The above consideration pertained to the limiting case c = 0, when the exact solutions 
were available. The next necessary step is to consider c 7^ 0, when no exact solution for the 
zero-velocity solitons could be found. We therefore started by using the known relaxation 
algorithm |jl5| in order, first of all, to obtain stationary soliton solutions numerically from 
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the ordinary differential equations produced by the substitution of Eqs. (j^) and (|^) into 
Eqs. (|1]) - (^. By properly setting boundary conditions, it was always possible to obtain 
a soliton solution for a given uj. 

A major objective here is to find out whether at fixed values of all parameters except c, 
there exists a maximum value of c above which the solitons are unstable. It has been found 
that, depending on the value of tu, there indeed exists Cmax beyond which solitons become 
unstable. However the instability at c > Cmax leads not to disappearance of solitons, but 
rather to their self-rearrangement into a slightly different form. 

A typical result is displayed in Fig. 2, with K=l,A = /i = and uo = 1.6. This 
value of UJ was chosen since it lies sufficiently deep inside the stability region at c = 0. 
Our simulations show that the value Cmax is very large, ~ 4.2. As is seen in Fig. 2b, at 
c > Cmax the soliton becomes unstable and, after shedding some radiation, it evolves into 
another member of the soliton family. 

A practically significant consequence of the above result is that values of c close to 1 
(recall that 1 group velocity in the nonlinear core) definitely give rise to stable solitons. 
This inference is important for experiments, because, in the most realistic case when both 
cores are made of the same material, the group velocities in them are necessarily close. 

In the case A 7^ (when BG is written on the linear core too), the relaxation algorithm 
also successfully generated stationary solitons. Starting with these, we have found that 
the solitons are stable in a broad parametric region, again including values of c essentially 
exceeding 1. However, detailed analysis of the interplay of A with other parameters is 
very cumbersome and is left aside. 

INTERACTIONS BETWEEN SOLITONS AND GENERATION OF MOVING 

SOLITONS 

Since the present model is nonintegrable, interactions between solitons may be quite 
complex. The simplest approach to simulating these interactions is to start from a su- 
perposition of two identical exact solitons placed initially at a distance from each other 
with some phase difference Aip. Results of the simulations, typical examples of which are 
displayed in Fig. 3, are similar for different values of the soliton's internal parameters and 
initial separation (provided that the solitons overlap weakly), but they strongly depend 
on Aip. In the case Aip = vr, the solitons, quite naturally, repel each other, cf. the well- 



known fact that nonlinear Schrodinger (NLS) solitons interact repulsively when A(y9 = vr 
||16|| . Even if the initial separation between the solitons is relatively large, the repulsion 
is strong enough to lend the two initially quiescent solitons conspicuous velocities, see 
Fig. 3a. In this case, eventual velocities are found to be W± = ±0.03. Thus, these 
simulations not only shed light on the character of the interaction between the solitons, 
but also provide a convenient way to generate stable moving ones. 

It is also interesting to compare the initial energy Ei of each soliton, defined as 
/-c» (l'"P + 1"^^ + I0P+ IV'P) d^ (which is a dynamical invariant of the model), and fi- 
nal values Ef of the energy of the moving solitons. In the case shown in Fig. 3a, 
Ef/Ei = 0.986, i.e., about 1.5% of the initial energy is lost (into emission of radiation) 
as a result of the interaction process. It should be stressed that moving solitons pro- 
duced by the interaction exhibit some internal vibrations, i.e., the solitons appear with a 
weakly excited internal mode (the existence of internal modes in stable GMTM solitons 



is a known fact |TT], |T2| ) . It may also happen that they capture some radiation which will 
be very slowly radiated away in the course of very long evolution (which is not relevant 
for experiments). 

In the opposite case Aip = 0, the solitons attract each other, which is similar to what is 
known for the NLS solitons. As is shown in Fig. 3b, they temporally merge into a single 
pulse, which later splits into two moving solitons with small internal vibrations. In this 
case, a conspicuous breaking of the initial symmetry between the two solitons is observed 
(special care has been taken to check that it is not an artifact produced by the numerical 
scheme). A plausible explanation is that the lump produced by the strong temporary 
overlapping of the initially attracting solitons (see Fig. 3b) is unstable against symmetry- 
breaking perturbations, the breaking being incomplete since the solitons separate quickly 
enough. This conjecture seems natural, as it is well known that various multisoliton states 



in the NLS equation are strongly unstable in the case of attraction |T^, but, of course, 
much more extensive simulations are necessary to check it in detail. 

This partial symmetry breaking can be characterized by the final/initial energy ratios 
for the two solitons shown in Fig. 3b, which is found to be Ef/E^ = 0.892 and 0.864 for 
the left and right solitons, respectively. In this case, a considerable share of the initial 
energy, ^ 12%, is lost into radiation. The final velocities of the solitons are W± = ±0.22, 
i.e., there is no tangible symmetry breaking in terms of the velocities. Note that \W±\ are 
much larger in this case than in the case Aip = it. 
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In the case Aip = 0, there is another noteworthy aspect of the symmetry breaking: 
the final sohtons demonstrate an internal asymmetry, characterized by the ratios of their 
partial energies, 

£_ = / |m_| dx/ / |t>_| dx, £+ = / |f+| dx/ / |m+| dx, 

J — OO J —CO J —CO J — CO 

where the positive (negative) subscript pertains to the right (left) soliton. For the solitons 
shown in Fig. 3b, this ratio takes values £_ = 0.470 and e+ = 0.465. In accord with these 
values, the final solitons, being intrinsically asymmetric, are, to a good approximation, 
mirror images of each other. 

We have also simulated the interaction of solitons with the initial phase difference is 
A(f = 7r/2. In this case the solitons repel each other, about 0.7% of the energy is lost into 
radiation, and the symmetry breaking is much more conspicuous, with the final velocities 
being W- = —0.023 and W.^. = 0.019. The stronger symmetry breaking in this case 
can be easily understood, as the symmetry of the initial configuration, which was taken 
as Usoi{x — ^xq) + iUso\{x + ^xq), xo being the initial separation between the solitons, is 
not compatible with the Eqs. (|l])-(^ and is therefore broken upon the propagation in a 
straightforward way. 

CONCLUSION 

In this paper, we have introduced a model consisting of two linearly coupled cores, one 
having the Kerr nonlinearity and the other being linear. A Bragg grating is written on the 
nonlinear core, while the linear one may or may not be equipped with the grating. The 
model allows to considerably extend the previously studied family of the Bragg-grating 
solitons. Exact solutions were found for zero-velocity solitons in a limiting case when the 
group-velocity terms are absent in the equations for the linear core, while in the general 
case solitons were found numerically. The main issue is their stability. We have found a 
nontrivial stability limit for them in terms of an internal parameter of the soliton family. 
Depending on the frequency u, the solitons may remain stable up to quite large values 
of the group velocity in the linear core. This strongly suggests that stable solitons can 
indeed be generated experimentally in dual-core systems, with the cores made of the same 
material. The vast stability region for the zero-velocity solitons in the dual-core model, 
found in this work, suggests a possibility to experimentally look for the corresponding 
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localized states with the fully trapped light. Interactions of initially separated solitons 
were investigated too, showing a considerable spontaneous symmetry breaking in the case 
when the solitons attract each other, which may be a result of a natural instability against 
symmetry-breaking perturbations. The interaction always results in appearance of stable 
moving solitons. 
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Figure Captions 

Fig. 1. Evolution of an asymmetrically perturbed soliton when c = (only the u 
-component is shown). The other parameters are X = fi = 0, k = 1, and the soliton's 
internal parameter (see Eqs. ^ is ^ = 7r/3. 

Fig. 2. Evolution of solitons at values of c slightly below and above Cmax^ (a) c = 4.1; 
(b) c = 4.3. The other parameters are n = 1, X = fi = 0, and uj = 1.6. 

Fig. 3. Interaction of two identical solitons with 6 = 7r/3, placed initially at a distance 
8 with two different values of the initial phase difference between the solitons: (a) Aip = it 
(repulsion); (b) Aip = (attraction). The other parameters are n = 1, X = fi = 0, and 
the evolution time is 400. 
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